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ABSTRACT 

Compact object mergers are one of the currently favored models for the origin of 
GRBs. The discovery of optical afterglows and identification of the nearest, presumably 
host, galaxies allows the analysis of the distribution of burst sites with respect to 
these galaxies. Using a model of stellar binary evolution we synthesize a population of 
compact binary systems which merge within the Hubble time. We include the kicks in 
the supernovae explosions and calculate orbits of these binaries in galactic gravitational 
potentials. We present the resulting distribution of merger sites and discuss the results 
in the framework of the observed GRB afterglows. 
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1 INTRODUCTION 



The recent discoveries ot A-ray a fterglows ot gamm a-ray 
bursts by the Beppo SAX satellite (Costa et al., 1997) have 



revolutionized the approach to these phenome na. For the 
first time s ince their discovery 30 years ago (Klebesadel 



et al., 1968) gamma ray bursts have been identified with 



sources at other wavelength s. In consequence o ptical after- 
glows have been discovered ( Groot et al., 1997a ), which lead 
to identification of host galaxies ( Groot et al., 1997b ). At 
the time of writing more than a dozen afterglows have been 
identified. The optical lightcurves of the GRB afteglows de- 
cay as a power law oc t~ a with the typical values of the 
index a = 1.1 to 1.3. In some cases the host galaxies have 
been found by observing the flattening of of the lightcurve. 
The underlying steady flux is identified as the emission of 
the host galaxy. In a few cases the host galaxies themselves 
were found as extended objects. This allowed to measure the 
offset between the location of the gamma-ray burst and the 
center of the host galaxy. A list five GRBs and the offsets 
from their host galaxies is shown in Table |l|. The offsets are 
generally small and the afterglows lie directly on the host 
galaxies. In other cases where the host galaxy have been 
only found from the flattening of the lightcurve we know 
that the location of the host galaxy and the GRB do not 
differ substantially from the simple fact that the host galax- 
ies have been identified. 

There are two basic categories of theoretical models of 
the central engines of gamma-ray bursts within the cosmo- 
logical model. The first class connects gamma-ray bursts 
with mergers of compact objects, e.g. neutron stars and/or 
black holes. There exist numerous scenarios in this set of 
models, some of them link GRB s with the coalescenc e of 
a black hole neutron star binary (Narayan et al., 1992). In 



other models l ike in the recent paper by (Kluzniak and Ru 



derman, 199£) the GRB events are related to mergers of 



two neutron stars. Compact object merger model provides 
enough energy to power a GRB, and it has been showed 



in the numerical simulations (Kluzniak and Lee, 1998) that 
the coalesce nce may last up to a se cond. The analytical esti- 
mations of (Portegies Zwart, 199£) extend this timescale up 
to a minute. The second class of models ( Paczynski, 1998] ) 



relates GRBs to explosions of supermassive stars, so called 
hypernovae. A direct prediction in this class of models is that 
gamma-ray bursts are related to the star forming regions. 

The relation between the GRB location and the host 
galaxies does not have to be true in the models that relate 
GRBs to compact object mergers. Tutukov and Yungelson 
( 1993[ ) have calculated the compact object merger rates, and 
also found that compact object binaries may travel the dis- 
tances up to 1000 kpc bef ore merging. In a more detailed 
study ( Bloom et al., 199*8^ ) calculated a population of com- 
pact object binaries using t he population synthesis method 
of Pols and Marinus ( 1994 ) and then calculated the spatial 



distribution of mergers in the potentials of galaxies for w 
few representative masses. They found that approximately 
15% of mergers take place outside the host galaxies). They 
have used the a Maxwellian kick velocity dist ribution with 
a v = 190 km s" 1 (Hansen and Phinney, 1997) 

In thi s work we use the the po pulation synthesis code 
based on ( Bethe and Brown, 1998] ) and extended by Bel- 
czynski and Bulik (199£). We concentrate on the depen- 



dence of the properties of the compact object binaries on 
the parameters used in the population synthesis code. We 
find that the most important parameter that determines the 
population of compact object binaries is the kick velocity a 
neutron star receives at birth, h oweve r this distribution is 
poorly known. Iben and Tutukov ( |l996| ) claim that the prop- 
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erties of pulsars can be explained by only the recoil veloci- 
t ies w ith no need for the kicks. Blaauw and Ramachandran 
(1998) find that a single kick velocity of 200 km s" 1 suffices 
t o rep roduce the pulsar population. Cordes and Chernoff 
(1997) proposed a weighted sum two Gaussians: 80 percent 
with the width 175 km s _1 and 20 percent with the width 
700 km s" 1 . 

We outline the model of the binary evolution and prop- 
agation in a galactic potential in section 2. The results of the 
calculation are presented in section 3 and we discuss them 
in section 4. 



2 MODEL 



2.1 Binary evolution 



In order to study the spatial distribution of compact object 
mergers we use the populatio n synthesis method. We use 
the population synthesis code (Belczynski and Bulik, 199!:) 
which concentrates on the population of massive star bi- 
naries, i.e. those that may eventually lead to formation of 
compact objects and compact object binaries. We include 
the evolution of the binaries due to interaction and mass 
transfer and also the kicks that a newly born neutron star 
receives in supernova explosion. A binary may be disrupted 
in each of the supernova events. The surviving binaries ob- 
tain center of mass velocities, which change their trajectories 
and may even eject them from their galaxy. 

While the evolution of single stars depends only on their 
mass and metallicity the evolution of binaries is also a func- 
tion the initial orbit (semimajor axis a, and eccentricity e) of 
the two stars. We assume that the distribution of the initial 
parameters can be expressed as a product of distributions 
of four parameters: the larger star (primary) mass M, the 
mass ratio of the less massive to the more massive star in 
the binary q, and the orbital parameters a and e, i.e that 
this quantities are in dependent. The distribut ion of primary 
masses used here is (Bethe and Brown, 1998) 

*(M) oc M" 3/2 , 

and we adopt a flat distribution of the mass ratio q. The 
semi major axis distribution is scale invariant, i.e. 

r(a) oc a -1 

with the limits 6Rq < a < 6000i?o , and we draw the eccen- 
tricity from a distribution H(e) = 2e. 

We assume that the kick velocity distribution is a three 
dimensional Gaussian, and parameterize it with its width 
cr„, i.e. 



f \ "±-3 2 u 

p(vj = —=o v v exp j 



(1) 



We generate population of compact object binaries for a few 
values of a v in order to asses the sensitivity of our results to 
this parameter. 

We describe the mass transfer in the common envelope 
evolution by the comm on envelope parameter acE (see e.g. 
(Vrancken et al., 1991)), and we use an intermediate value 



of 0.8 for this parameter. In this type of evolution the more 
massive star looses its envelope and becomes a helium star 



Table 1. GRBs with measured offsets from the centers of their 
host galaxies. 

Reference 



GRB 


redshift z 


Offset A0 


970228 


??? 


0.30" 


970508 


0.835 


0.01" 


971214 


3.42 


0.06" 


980703 


0.966 


0.21" 


990123 


1.60 


0.60" 



(Sahu et al.. 1997) 



f Fruchter. 1998j1 

( Kulkarni et al.. 199g| ) 



(Bloom et al.. 1998b) 



(Djorgovski, 1999) 



with mass approximately 30% of its initial value. The (3 pa- 
rameter which describes the specific angular momentum of 
the material expelled from the binary in the Roche lo be over- 
flow phase is set to ft = 6 (Pols and Marinus, 1994). Accre- 
tion onto a neutron star in a binary is treated as B ondi-Hoyle 



accretion and we use the formalism developed by ( Bethe and 



Brown, 1998) to find the amount of mass accreted onto the 
neutron star, and the final orbital separation. Systems with 
nearly equal masses evolve at the similar speed, and loose 
the common envelope, shrinking their orbit at the same time. 
For a mo re detailed description of th e population synthesis 
code see (Belczynski and Bulik, 1999). 

We assume that a neutron star with mass of 1.4 Mq 
is formed in each supernova explosion. We draw a random 
time in the orbital motion to obtain the position on the orbit 
when the supernova explodes. The remaining mass of the 
envelope is ejected from the system, and the newly formed 
neutron star receives a kick. We verify whether the system is 
still bound after the explosion. For bound systems we find 
the parameters of the new orbit and the kick velocity the 
whole binary receives 



AV 



Mi 



Mi + M. 



Ml 

J (V 2 + Vkick ) 



Ml are the 



where M\ is the mass of the companion, M\ 
initial and final masses of the supernova, vi is the orbital 
velocity of the supernova at the time of explosion. After each 
supernova expolsion we verify whther the system survives as 
a binary. 

A compact object binary loses its ene rgy through g rav- 
itational radiation. The time to merge is (Peters, 1964) 



5cV(l - 



256G s Mm(M + ; 



' , 73 2 37 4 
, 1+ 24 C + 96 6 



(2) 



where a is the semi major axis of the orbit, e is its eccen- 
tricity, and M, m are the masses of the compact objects. 



2.2 Orbit in Potentials of Galaxies 

Since little is known about the host galaxies of gamma- 
ray bursts, in particular of their types and masses, we will 
present two extreme cases: (i) propagation in a potential of 
large spiral galaxy like the Milky Way, and (ii) propagation 
in empty space, corresponding to GRBs originating e.g. in 
globular clusters. In the latter case we assume that all bina- 
ries originate in one point, and travel due to kicks described 
above and there is no gravitational potential. 

The potential of a spiral galaxy can be described as 
a sum of three components: bulge, disk, and dark matter 
halo. A convenient wa y to describe the Galactic potential 
has been proposed by ( Miyamoto and Nagai, 1975 ), while a 
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Figure 1. Distribution of compact object system in the velocity versus lifetime (time to merge). The top left panel shows the case when 
there is no kick velocity, the top right panel shows the case cr = 200km s — 1 , the bottom left panel C7 = 400km s -1 , and the bottom right 
panel is for <r = 800km s~ x . The horizontal dashed line corresponds to the Hubble time (15Myrs). In the region for t merge < 15M we 
present two solid lines: the vertical corresponding to v = 200 km s _1 - approximately the escape velocity from a galaxy, and the line 
corrsponding to a constabt value of v X tm erge = 30 kpc. Together these lines define the region in the parameter space with systems that 
can escape from the host galaxy. 



and Gi 



series o f more detai led models were constructed by (Kuijkcn 
more, 1989| ) and used in modeling; the galact i c halo 



popula tion of neu tron stars (Bulik and Lamb, 1995 ; Bulik 
et al., 1^98). The ( |Miyamoto and Nagai, 1975| ) potential for 
a galactic disk and bulge is 



GM 



where a,j and b j are the parameters, M is the mass, and 
R = y/ x 2 = y 2 . The dark matter halo potential is spheri- 
cally sy: 

$( r ) = 



cally symmetric 
GM h 



corresponds to a mass distribution p — p c /[l + (r/r c ) 2 ]. 
The mass of such halo is infinite, so we introduce a cutoff 
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Figure 2. The cumulative distributions of the projected dis- 
tances ri of compact object binaries in a potential of a large 
galaxy. The curves correspond to the four values of the kick 
velocityo- = 0, 200, 400, 800km s — 1 . The vertical line corrsponds 
to the distance of 30 kpc from the center of galaxy. 



0.8 - 



0.4 



0.2 



1 1 1 

kms- 1 


1 1 


200 kms- 1 


/--•/ 


400 kms" 1 

800 kms" 1 


: / 


/// ; 












1 1 1 





10" 1 



10° 



10> 10 2 
rjkpc] 



10 3 



10 4 



Figure 3. The cumulative distributions of the projected 
distances rj_ of compact object binaries in empty space. 
The curves correspond to the four values of the kick 
velocitycr=0, 200, 400, 800 km s .The vertical line corrsponds to 
the distance of 30 kpc from the point where the systems originate. 



value r cu t = 100 kpc above which the density of the halo 
falls to zero. While the details of the model of galactic po- 
tential are not important for this study we have to adopt 
a particular value of the masses and sizes of each of the 
component s. We use the values of the parameters as deter- 
mined by ( Blacs and Rajagopal, 1991 ) for the Milky Way: 
oi = Okpc, bi = 0.277 kpc, a 2 = 4.2 kpc, b 2 = 0.198 kpc, 



Mi = 1.12 x 1O 1U M , M 2 = 8.78 x 1O 1U M , r c = 6. Okpc, 
and M h = 5.0 x 10 10 M Q . 

We assume that the distribution of binaries in our 
model galaxy foll ows the mass distribution in the young disk 
(Paczyriski, 199C), that is 



P(R, z) dRdz = P{R)dRp{z)d(z) . 

The radial distribution is exponential with 

P(R) oc e- R/Re **R 



with R e 



■■ 4.5 kpc and extends to R„ 



tical distribution is p(z) oc e 



-z/z e 



and 



20 kpc. The ver- 
ex P = 75 pc. We 



note that this is not a self consistent approach: the density 
inferred from the disk potential is not the same as the den- 
sity of binaries. However in this work we are not interested in 
determining high accuracy positions around the host galaxy, 
and rather with an estimate of the general properties of the 
distribution of compact object mergers. 

Each binary moves initially with the local rotational 
velocity in the galactic disk. After each supernova explosion 
we add an appropriate velocity, provided that the system 
survives the explosion. We calculate the orbit of each system 
until it merger time provided that the merger time is smaller 
than the Hubble time (15 Gyrs here). 



3 RESULTS 

The kick velocity distribution is not very well known. There- 
fore, we use the population synthesis code with four values 
of the kick velocity distribution width: with no kick veloc- 
ities a v = km s" 1 , and with a v = 200, 400, 800 km s" 1 . 
This covers the range of values this distribution is likely to 
have. This the same approach as adopted in our previous 
work ( Belczynski and Bulik, 199£ ). 

The binaries receive kicks for two reasons. First, the 
envelope of the supernova is lost from the system and it 
carries away some momentum. Thus even in the case when 
there is no kick veloci ty a binary achieves an additional ve- 
locity (Blaauw, 196C). Second, if the supernova explosion 
is asymmetric both the newly formed compact object may 
receive a kick velocity which affects the orbit of the binary 
after the explosion as well as its center of mass velocity. The 
fate of a binary system in a supernova explosion depends on 
the value and direction of the kick velocity, on the orbital 
phase at which the explosion occurs, and on the parameters 
of the binary: the masses and orbital parameters a, and e. 

We present the population of compact object binaries 
in the plane spanned by the center of mass velocity ofter the 
second supe rnova expolsio n and time to merge in Figure |l| 
The orbital (Blaauw, 1960) effects are isolated and shown in 
the top left panel of Figure [j], where we present the results 
of the simulation with a v — 0. There is a tail of long lived 
systems with lifetimes much longer than the Hubble time 
and small velocities, stretching outside of the boundaries of 
the plot to the lifetimes even of 10 20 years. These systems 
originally had large orbital separations, and hardly inter- 
acted in the course of their binary lifetime. In the case when 
there are no kicks the center of mass velocity of the comapct 
object binary depends on the amount of mass lost in the su- 
pernova explosion. In the extreme case of large mass loss, 
the center of mass velocity approaches the orbital velocity 
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Figure 4. The distribution of compact object mergers around a small galaxy (approximated by empty space). The region shown in the 
plots increases from 10 kpc in the top left panel, through 100 kpc in the top right panel, to 1000 kpc in the bottom left panel and to 
finally to 10 Mpc in the bottom right panel. 



at the moment of supernova explosion, and it can never ex- 
ceed it. The velocity of the system increases with increasing 
mass loss, however the systems that loose too much mass 
become unbound. This is why the lower part of the plot be- 
low t m erge ~ 10 8 years is empty. With increasing the kick 
velocity also the typical velocity of a system increases and 
there appear short lived systems in tight orbits. They can 
now survive a large mass loss when the kick velocity has a 
favorable direction. Thus as the kick velocity is increased 
only the tightly bound systems (with short merger time) 
can survive the supernova explosions. Another effect of the 
kick velocity is that the long lived systems with t merge much 
longer than the Hubble time, which were present in the case 
a v — km s _1 disappear. The typical velocity of a system in- 
creases with the kick velocity. However, the population of the 



comapct merger binaries is not much affected when the kick 
velocity becomes large, e.g. changing the kick velocity dis- 
tribution width from a v = km s _1 to 200 km s _1 produces 
a much stronger effect than going from a v = 400 km s _1 to 
800 km s _1 . Most of the systems are disrupted by such high 
velocities, and the surviving ones are only those for which 
the kick are not so large and have a favorable direction. 

Another effect of increasing the kick velocity is that 
the typical lifetime of a system becomes smaller. When the 
kick velocity is large only very tight, and/or highly eccen- 
tric systems survive, hence the typical lifetime of compact 
object binaries decreases. It should be noted the typical cen- 
ter of mass velocity of the compact object binaries increases 
roughly linearly with the kick velocity, while the lifetime 
decreases approximately exponentially. 
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Figure 5. The distribution of compact object mergers around a massive galaxy. The region shown in the plots increases from 10 kpe in 
the top left panel, through 100 kpe in the top right panel, to 1000 kpe in the bottom left panel and to finally to lOMpc in the bottom 
right panel. 



In Figure |l] we also plot following (Bloom et al., 1998a) 
the lines corresponding to the Hubble time (the dashed line), 
and we mark the region with the stars that will escape from 
a galactic potential. In order to escape a binary must satisfy 
the following conditions: (i) it has to have a velocity larger 
than the escape velocity, (ii) the distance wt mer9e must be 
larger than the size of the galzxy. We also draw the line at 
tmrg = 15Myrs, to denote the systems that merger within 
the Hubble time. All the systems to the right of the solid 
line in Figure refvt have velocities above 200 km s , and 
live long enough to travel further than 30 kpe. We should 
also note that although each panel in Figure |l| contans 10 3 
systems, the production rate of compact object binaries de- 
creases exponentially with the i ncrea sing kick velocity (see 
eq. 13 in Belczyriski and Bulik (|l999[)). 



In Figures |^ and |§| we present the cumulative distribu- 
tions of the projected distance from the center of the host 
galaxy in case (i) and from their the place of birth in case 
(ii), respectively, of the systems that merge within the Hub- 
ble time. When the binaries propagate in the potential of 
a large galaxy the kick velocity only weakly influences the 
the distribution of the mergers. Below the radius of 10 kpe 
the distribution is determined by the potential well. This 
is where all short lived and slow systems merge. There ex- 
ists however a tail of high velocity, long lived systems (see 
Figure that manage to escape. The escaping fraction is a 
weak function of the kick velocity. Typically the number of 
systems that merge further than 30 kpe from the center of 
the host galaxy is 30%, except for the unphysical case of no 
kick velocities when it drops below 20%. 
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In the other extreme case of small host galaxies for 
which we neglect the gravitational potential the escaping 
fraction can be even larger. The escaping fraction decreases 
from 80% for the kick velocity a v = 200 km s" 1 to about 
50% for a v = 800 km s" 1 . The reason for such behavior is 
clearly seen from Figure |lj In the product of center of mass 
velocity and time to merge the dominant role is played by 
the fast decrease of the time to merge with the increasing 
kick velocity a v . 

These quantitative results are visualized in Figures ^ 
and ^. Here we show the distribution of 10 3 mergers around 
massive galaxy and in the empty space. We are showing four 
panels that cover the scales from 10 kpc to 10 Mpc. In the 
case of the propagation in a massive galaxy we are showing 
the projection in the plane of the galactic disk so the effects 
of the r otational velocity and the asymmetry of the potential 



well are visible. Both calculations have been done for the 



case of the kick velocity distribution width a v = 200 km s 



4 CONCLUSIONS 

We find that a significant fraction i.e. more than 20 percent 
of compact object mergers take place outside of the host 
galaxies. The figures obtained for our case of no gravita- 
tional potential should be considered as an upper limit only. 
In contrast the observatio ns show that the GRB aft erglows 
lie on the host galaxies ( Hogg and Fruchter, 1998 ). How- 
ever, our sample of the observed GRBs with afterglows is 
yet limited to the bursts longer than 6 s as BeppoS AX trig- 
gers on this timescale. Long bursts could be connected with 
the hypernovae-like events and therefore they are closely as- 
sociate d with the galaxies . Compact object mergers may be 



con needed with the short hursts although it. has heen argued 
that mass transfer in the coalescence of c ompact object may 
last much longer (Portegies Zwart, 1998). 

Our results are consistent with the calculation by 



(Bloom et al., 1998a), for the case of a massive galaxy and 
the kick velocity distribution width a v = 200 km s~ . We 
have verified the dependence of the distribution of compact 
object mergers on a v . Our results show that for the case 
of massive galaxies the escaping fraction weakly depends 
on the distribution of kick velocities. We include also bina- 
ries with objects that are higher mass than the canonical 
1.4 M©. These binaries are formed through accretion from a 
giant companion onto the neutron star. In this calculation 
the highest mass of a compact object is below 2.5 Mq. The 
distribution of these more massive binaries is slightly more 
concentrated around the galaxies. 

It has to be noted that there is a numnber of poten- 
tial selection effects which may affect the results of this 
study. Assuming that compact objcet mergers are reponsible 
for GRBs there may be qualitative differences between the 
NS-NS mergers and NS-BH mergers. As indicated above, 
their spatial distribution around host galaxies is different. 
Also the typical timescale of the bursts may be different 
between these two classes. Gamma-ray bursts form two sep- 
arate classes (long vs. short) with different brightness dis- 
tributions and spectra (short burst are harder than the long 
ones). So far the study of afterglows has been possible only 
for the long bursts. It may be the case that compact object 



mergers are connected with the short bursts for which so far 
no information about the host galaxies exist. 

The host galaxies have been identified in a long obser- 
vational procedure: a gamma-ray burst lead to identification 
of a fading X-ray source, and then to discovery of the opti- 
cal afterglow. Precise observations of the optical afterglows 
lead to the disory of host galaxies. There are bursts for which 
the X-ray or optical afterglows were not found. Since after- 
glows are usually connected with external shocks, gamma- 
ray bursts that take place outside of galaxies have much 
weaker afterglows because of the low density of int ergala ctic 
matter. Begelmann etal. Begelman and Meszaros (1993) ar- 
gue that the afterglow emission depends scales only with the 
square root of the density of the outside medium. The mean 
external densities measured from the analysis of th e known 
afterglow lightcu rves are typically n ~ 0.03 cm~ 3 ( ]Galama 



and Wijers, 1998), while the intergalactic medium may be as 
rarificd as 10 _t ' cm -3 . Hence the afterglow of a burst taking 
place outside a galaxy may be up to two orders of magnitude 
weaker than the one in a galaxy. It shows that there may 
be a strong preference against identification of host galaxies 
for the bursts that take place outside of galaxies. 
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